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The singular dynamics at critical attractors of even the simplest one di- 
mensional nonlinear iterated maps is of current interest to statistical physicists 
because it provides insights into the limits of validity of the Boltzmann-Gibbs 
(BG) statistical mechanics. This dynamics also helps inspect the form of the 
possible generalizations of the canonical formalism when its crucial supports, 
phase space mixing and ergodicity, break down. 

The fame of the critical attractors present at the onset of chaos in the logistic 
and circle maps stems from their universal properties, comparable to those of 
critical phenomena in systems with many degrees of freedom. At these attractors 
the indicators of chaos withdraw, such as the fast rate of separation of initially 
close by trajectories. As it is generally understood, the standard exponential 
divergence of trajectories in chaotic attractors suggests a mechanism to justify 
the property of irreversibility in the BG statistical mechanics P^. In contrast, 
the onset of chaos imprints memory preserving properties to its trajectories. 

The dynamical nature of trajectories is appraised on a regular basis through 
the sensitivity to initial conditions £t, defined as 

ft(»o)= lim (Ax t /Axo), K t, (1) 

Ai - >0 

where Axq is the initial separation of two trajectories and Axt that at time t. 
For a one-dimensional map it has the form £t(cco) = exp(Aii), with Ai > for 
chaotic attractors and Ai < for periodic ones. The number Ai is called the 
Lyapunov coefficient. At critical attractors Ai = and £t does not settle onto 
a single- valued function but exhibits instead fluctuations that grow indefinitely. 
For initial positions on the attractor £ t develops a universal self-similar temporal 
structure and its envelope grows with t as a power law. 

It has been recently corroborated 0-0 that the dynamics at the critical 
attractors associated to the three familiar routes to chaos, intermittency, period 
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doubling and quasiperiodicity 0, obey the features of the g-statistics, the gen- 
eralization of BG statistics based on the g-entropy S q 0. The focal point of the 
g-statistical description for the dynamics at such attractors is a £t associated to 
one or several expressions of the form 

£t(xo) = exp q [X q (x ) t], (2) 

where q is the entropic index and X q is the g-generalized Lyapunov coefficient. 
Also the identity K\ — X\ (where the rate of entropy production K\ is given by 
K±t = Ssait) ~ Sbg(0) with Sbg the Boltzmann-Gibbs entropy) generalizes to 

K q = X q , (3) 

where the rate of g-entropy production K q is defined via K q t — S q (t) — S q (0) 
0,0,111. 



Tsallis q index & Mori's g-phase transitions 

The central issue of research in g-statistics is perhaps to confirm the oc- 
currence of special values for the entropic index q for any given system and to 
establish their origin. In the case of critical attractors the allowed values for q 
are obtained from the universality class parameters to which the attractor be- 
longs. For the simpler cases, the pitchfork and the tangent bifurcations, there 
is a single well-defined value for the index q 0. The pitchfork bifurcations 
form the sequence of period doublings that culminate in the chaos threshold 
whereas at the tangent bifurcation chaos develops via intermittency. For crit- 
ical multifractal attractors, as in the period doubling accumulation point and 
in the quasiperiodic onset of chaos, the situation is more complicated and there 
appear to be a multiplicity of indexes q but with precise values given by the 
attractor scaling functions. They come out in pairs and are related to the oc- 
currence of dynamical 'g-phase' transitions 0, and these are identified as the 
source of the special values for the entropic index q. The g-phase transitions 
connect qualitatively different regions of the multifractal attractor. 

The main quantities in the thermodynamic formalism of g-phase transitions 
(developed by Mori and colleagues in the late 80's 0) are the spectral functions 
<p(q) and ipW, related to each other via Legendre transformation, the func- 
tion of generalized Lyapunov coefficients A(q), given by A(q) = d0(q)/dq, and 
the variance v(q) = dX(q)/dq. The functions <^(q) and ^(A) are the dynamic 
counterparts of the multifractal dimensions -D(q) and the spectrum f(a) that 
characterize the geometric structure of the attractor jSj . As with ordinary ther- 
mal 1st order phase transitions, a g-phase transition is indicated by a section of 
linear slope m c = 1 — q in the 'free energy' ip(\), a discontinuity at q = q in the 
'order parameter' A(q), and a divergence at q in the 'susceptibility' v(q). Actu- 
ally an infinite family of such transitions takes place but of rapidly decreasing 
strengths 0, 0. 
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• Tangent & pitchfork bifurcations 

The tangent bifurcations of unimodal (one hump) maps of general nonlin- 
earity z > 1 display weak insensitivity to initial conditions, i.e. power-law 
convergence of orbits when at the left-hand side (x < x c ) of the point of tan- 
gency x c . However at the right-hand side [x > x c ) of the bifurcation there is a 
'super-strong' sensitivity to initial conditions, i.e. a sensitivity that grows faster 
than exponential The two different behaviors can be couched as a q-phase 
transition with indexes q and 2 — q for the two sides of the tangency point. The 
pitchfork bifurcations display weak insensitivity to initial conditions. 

For the tangent bifurcations one has always q = 3/2, while for the pitchfork 
bifurcation one has q = 5/3. Notably, these results for the index q are valid for 
all z > 1 and therefore define the existence of only two universality classes for 
unimodal maps 2 . The treatment of the tangent bifurcation differs from other 
studies of intcrmittency transitions in that there is no feed back mechanism 
of iterates from an adjacent chaotic region. Therefore, impeded or incomplete 
mixing in phase space (a small interval neighborhood around x — x c ) arises 
from the special 'tangency' shape of the map at the transitions that produces 
monotonic trajectories. This has the effect of confining or expelling trajectories 
causing anomalous phase-space sampling, in contrast to the thorough coverage 
in generic states with Ai > 0. 

• Period-doubling accumulation point 

For a unimodal map of nonlincarity z > 1 (e.g. the logistic map has z = 2) 
with extremum at x = and control parameter jj, the onset of chaos is obtained 
at the accumulation point \ioo of the \x values for the pitchfork bifurcations \i n , 
n = 1, 2, [H]. This is often called the Feigenbaum attractor which reappears 
in multiples together with the precursor cascade of period-doubling bifurcations 
in the infinite number of windows of periodic trajectories that interpose the 
chaotic attractors beyond /ioo. The number of cascades within each window 
is equal to the period of the orbit that emerges at the tangent bifurcation at 
its opening. See Fig. 1. The dynamics at the Feigenbaum attractor has been 
analyzed recently 0]- By taking as initial condition xq = at fioo it is found 
that the resulting orbit of period 2°° consists of trajectories made of intertwined 
power laws that asymptotically reproduce the entire period-doubling cascade 
that occurs for \i < /ioo. See Fig. 2a. It was established that £t has precisely 
the form of a set of interlaced q-exponentials, of which the q-indexes and the 
sets of associated q-Lyapunov coefficients X q were determined. As mentioned, 
the appearance of a specific value for the q index (and actually also that for 
its conjugate value Q = 2 — q) turns out to be due to the occurrence of q- 
phase transitions between 'local attractor structures' at fx^. The values of 
the q-indexes are obtained from the discontinuities of the universal trajectory 
scaling function a. This function characterizes the multifractal by measuring 
how adjacent positions of orbits of period 2 n approach each other as n —> oo 
PJ. The main discontinuity in a is related to the most crowded and most sparse 
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regions of the attractor and in this case q = 1 — ln2/(^ — 1) Inap, where qj? is 
the universal scaling constant associated to these two regions. Furthermore, it 
has also been shown 0] that the dynamical and entropic properties at /itoo 
are naturally linked through the ^-exponential and g-logarithmic expressions, 
respectively, for the sensitivity £ t and for the entropy S q in the rate of entropy 
production K q . 

• Quasiperiodicity & golden mean route to chaos 

A recent study [5] of the dynamics at the quasiperiodic onset of chaos in 
maps with zero slope inflection points of cubic nonlinearity (e.g. the critical 
circle map) shows strong parallelisms with the dynamics at the Feigenbaum 
attractor described above. Progress on detailed knowledge about the structure 
of the orbits within the golden-mean quasiperiodic attractor, see Fig. 2b, helped 
determine the sensitivity to initial conditions for sets of starting positions within 
this attractor jS] . It was found that £ t is made up of a hierarchy of families of 
infinitely many interconnected g-exponentials. Here again, each pair of regions 
in the multifractal attractor, that contain the starting and finishing positions 
of a set of trajectories, leads to a family of g-exponentials with a fixed value of 
the index q and an associated spectrum of g-Lyapunov coefficients \ q . 

As in the period doubling route to chaos, the quasiperiodic dynamics con- 
sists of an infinite family of g-phase transitions, each associated to trajectories 
that have common starting and finishing positions located at specific regions of 
the multifractal. The specific values of the variable q (in the thermodynamic 
formalism) at which the g-phase transitions take place are the same values for 
the entropic index g in £t. The transitions come in pairs at g and 2 — g as 
they are tied down to the expressions for \ q in £ t . Again, the dominant dy- 
namical transition is associated to the most crowded and sparse regions of the 
multifractal, and the value of its g-index is [5] q = 1 — \nw gm /2 lna sm where 
w gm = (Vo— l)/2 is the golden mean and a gm is the universal constant that 
plays the same scaling role as Qf. 

• Structure in dynamics 

The dynamical organization within critical multifractal attractors is difficult 
to resolve from the consideration of a straightforward time evolution, i.e. the 
record of positions at every time t for a trajectory started at an arbitrary position 
Xq within the attractor. In this case what is observed are strongly fluctuating 
quantities that persist in time with a scrambled pattern structure that exhibits 
memory retention. Unsystematic averages over xo would rub out the details 
of the multiscale dynamical properties we uncovered. On the other hand, if 
specific initial positions with known location within the multifractal are chosen, 
and subsequent positions are observed only at pre-selected times, when the 
trajectories visit another selected region, a distinct g-exponential expression for 
£t is obtained. 
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Manifestations of q-statistics in condensed matter problems 

There are connections between the properties of critical attractors referred-to 
here and those of systems with many degrees of freedom at extremal or transi- 
tional states. Three specific examples have been recently developed. See Table. 
In one case the dynamics at the chaos threshold via intermittency has been 
shown to be related to that of intermittent clusters at thermal critical states 
|1(J| . In the second case the dynamics at the noise-perturbed period-doubling 
onset of chaos has been found to show parallelisms with the glassy dynamics 
observed in supercooled molecular liquids In the third case the known con- 
nection between the quasiperiodic route to chaos and the localization transition 
for transport in incommensurate systems is analyzed from the perspective of 
the q-statistics \12\ . 

• Critical clusters & intermittency 

The dynamics of fluctuations of an equilibrium critical state in standard 
models of a magnetic or fluid system is seen to be related to the dynamics at 
a critical attractor for intermittency. The connection can be examined when 
instead of the entire critical system a single unstable cluster of excess magne- 
tization or density of size R is considered. The analysis, initially developed by 
Contoyiannis and colleagues |13| . has been reconsidered recently in connection 
to q-statistics [Tfj] . 

An important element in the analysis is the determination of the order pa- 
rameter 4>{r) of a large cluster by withholding only its most probable configu- 
rations from a coarse-grained partition function Z . The conditions under which 
these configurations dominate are evaluated as these determine an instability 
of the cluster. The instability can be expressed as an inequality between two 
lengths in space. This is ro R, where r$ is the location of a divergence in 
the expression for 4>{r). In a coarse-grained time scale the cluster is expected 
to evolve by increasing its average amplitude (j> and/or size R because the sub- 
system studied represents an environment with unevenness in the states of the 
microscopic degrees of freedom (e.g. more spins up than down). Increments in 
<p for fixed R takes the position ro for the singularity closer to R, the dominance 
of this configuration in Z decreases accordingly and rapidly so. When the diver- 
gence is reached at ro = R the profile (j>(r) no longer describes the spatial region 
where the subsystem is located. But a subsequent fluctuation would again be 
represented by a cluster 4>{r) of the same type. From this renewal process we 
obtain a picture of intermittency ^U] ■ 

Amongst the static and dynamical properties for such single critical cluster 
of order parameter 4>(r) we mention |10j : 1) The faster than exponential growth 
with cluster size R of the space-integrated 4> suggests nonextensivity of the BG 
entropy but extensivity of a g-entropy expression. 2) The finding that the time 
evolution of <p is described by the dynamics of the critical attractor for inter- 
mittency which implies an atypical sensitivity to initial conditions compatible 
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with g-statistics. 3) Both, the approach to criticality and the infinite-size cluster 
limit at criticality manifest through a crossover from canonical to g-statistics. 

• Glassy dynamics &: noise-perturbed Feigenbaum attractor 

The erratic motion of a Brownian particle is usually described by the Langevin 
theory. As it is well known, this method finds a way round the detailed con- 
sideration of many degrees of freedom by representing the effect of collisions 
with molecules in the fluid in which the particle moves by a noise source. The 
approach to thermal equilibrium is produced by random forces, and these are 
sufficient to determine dynamical correlations, diffusion, and a basic form for 
the fluctuation-dissipation theorem. In the same spirit, attractors of nonlinear 
low-dimensional maps under the effect of external noise can be used to model 
states in systems with many degrees of freedom. In a one-dimensional map with 
only one control parameter [i the consideration of external noise of amplitude 
a could be thought to represent the effect of many other systems coupled to it, 
like in the so-called coupled map lattices The general map formula can be 
seen to represent a discrete form for a Langevin equation with nonlinear 'friction 
force' term 

The dynamics of noise-perturbed logistic maps at the chaos threshold has 
been shown to exhibit parallels with the most prominent features of glassy dy- 
namics in, for example, supercooled liquids. In this analogy the noise amplitude 
a plays a role similar to the temperature difference from a glass transition tem- 
perature. Specifically our results are [TTj: 1) Two-step relaxation occurring in 
trajectories and in their two-time correlations when a — > 0. 2) A map equivalent 
to a relationship between the relaxation time and the configurational entropy 
3) Both, trajectories and their two-time correlations obey an 'aging' scaling 
property typical of glassy dynamics when a — 0. 4) A progression from normal 
diffusiveness to subdiffusive behavior and finally to a stop in the growth of the 
mean square displacement as demonstrated by the use of a repeated-cell map. 
See Fig. 3. The existence of this analogy is perhaps not accidental since the 
limit of vanishing noise amplitude a — > involves loss of ergodicity. 

• Localization & quasiperiodic onset of chaos 

One interesting problem in condensed matter physics that exhibits connec- 
tions with the quasiperiodic route to chaos is the localization transition for 
transport in incommensurate systems, where the discrete Schrodinger equation 
with a quasiperiodic potential translates into a nonlinear map known as the 
Harper map |14j . In this equivalence the divergence of the localization length 
corresponds to the vanishing of the ordinary Lyapunov coefficient. It is inter- 
esting to note that the basic features of g-statistics in the dynamics at critical 
attractors mentioned here turn up in the context of localization phenomena. 

Acknowledgements. Partially supported by CONACyT and DGAPA- 
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Route to chaos 


Intermittency 


Period doubling Quasiperiodicity 


Common 
properties 


Vanishing ordinary Lyapunov coefficient, 
dynamical phase transitions (Mori's g-phases), 
power-law dynamics, g-sensitivity, g-Pesin identity 


Distinctive 
properties 


(Also) faster 
than exponential 
dynamics 


Foam-like 
phase space 


Dense 
phase space 


Applications in 
condensed 
matter physics 


Critical 
clusters 


Glass 
formation 


Localization 


Applications in 

other 

disciplines 


Information & 
other flows 
in networks,. . . 


Protein folding, 
vegetation 
patterns,. . . 


Mode locking, 
cardiac cells, 
Internet TCP 



Table 1: Summary. The three routes to chaos, properties and applications. 
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Figure 1 : (a) The classic logistic map attractor as a function of control parameter 
fi. (b) Enlargement of the box in (a). 
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Figure 2: (a) Absolute values of positions in logarithmic scales of the first 10000 
iterations for the trajectory of the logistic map at the onset of chaos /ioo with 
initial condition xq = 0. The numbers correspond to iteration times. The power- 
law decay of time subsequences mentioned in the text can be clearly appreciated, 
(b) Positions 9(t) vs t in logarithmic scales for the orbit with initial condition 
9(0) = of the critical circle map for the golden-mean winding number. The 
labels indicate iteration time t where iffl is the Fibonacci number of order n. 




Figure 3: Glassy diffusion in the noise-perturbed logistic map. (a) Repeated- 
cell map (blue) and trajectory (red), (b) Time evolution of the mean square 
displacement (x^ \ — (x t ) for an ensemble of 1000 trajectories with initial con- 
ditions randomly distributed inside [—1,1]. Curves are labeled by the value of 
the noise amplitude. 
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